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Introduction 

Nitrous oxide is present in the atmosphere in minute quantities 

fnnhTT ' T l ° n u in 1988 WaS ~ 30? parts per billion b y volume 
(ppbv), about a thousand times less than that of COo and it is 

increasing at the rate of ~0.8 ppbv/yr (Elkins et al„ in press) The 
seemingly small growth rate, -0.25%/yr, is the result of a large 
imbalance (-30%) between the sources and sinks. S 

Despite its low abundance in the atmosphere, N 2 0 plays an impor- 
tant role in the radiative and chemical balance of the atmosphere. The 
e remey ong lifetime of N 2 0, -150 years, means the system has a 
veiy long memory of its emission histoiy. The radiative forcing of N 9 0 

uto) N O is r’r^H iS ; 2 °° UmeS ^ ° f C ° 2 ^ughton et ^ 
tinn° h 1 2 °, destr °y ed ln stratosphere by photolysis and by reac- 
tion with electronically excited oxygen atoms (0(i D )], making it the 

ominant precursor of odd nitrogen (NO x ) in the stratosphere. Thus 
increase in N 2 0 would lead to an increase in stratospheric NO 
which would catalyze the destruction of stratospheric ozone. 

.JPf f° Urces of N 2° are n °t well known. Of the -15 Tg N (1 Tg = 
10 kg) produced annually, by far the largest source seems to be 

fmm ,h nS natUra ' S °" S ’ 6 * 3 Tg/>T ’ f0ll0Wecl ^ missions 

strength of these natural sources is estimated from only a handful 

bustmn m h aSUrem T tS ' ° ther s ° urces lncIude emlsslon com- 
b, T“ S burning, and agricultural Helds amended with 
nitrogenous fertilizers. For some time, it was thought that, like CO„, 
p mary cause for the Increasing trend was combustion of fossil 
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fuels in particular, coal-burning power plants producing electacrty 
mio et al 19871. However, recent identification of an artifact in the 
S sampling procedure ruled out combustion (Including biomass 
burning) as the major cause of the N 2 G trend (Murto and KramUch. 
1988). We do not know why N 2 0 is increasing in t e a ™° s P ^ 
Global models of the N a O cycle (e.g., Prinn et al ’ 1990) 
tions in the concentrations of N 2 G in the atmosphere t^infe rfta 
location and magnitude of the N 2 G sources ^ ***^ ^_ 
point to large source(s) located in the tropics that ^"tr^ute : 

Lntlv to the latitudinal gradient and secular trend of N z O in the 
a^lhere SuTmi approach, however, cannot distinguish among 
individual sources without direct information about at least some o 

th ^n ^separate paper iBouwnan et al., in preparation), we focused, 

«* aSSEp budget itself- but on , **££££%££ 

in the budget i.e., emissions from natural soils. The eilects 
izer application and other human perturbations on the N 2 0 em- 
ails from soils were not considered. A very simplistic model was 
developed for estimating directly the geographic and seasonal varia 
S ^ 3 P f N 2 0 emissions from natural soils. We hope that the model, 
and the sensitivity studies using it, will provide guidance for mea- 
^femenrstrategTes to reduce systematically uncertainties abou 

^anllyses. Here, we present an account of the assumptions and 
decisions involved in the construction of the model. 


A Conceptual Model of N 2 0 Production 

N O is produced by both nitrification and denitrification. Demtri- 
flcatmn S occurs un/er oxygen limiting -ndlt.ons^. = us 

Davidson 1989), Under strictly anaerobic conditions nitrous oxide 

(1986) for a Review of the biogeochemistiy and measurement 

N 2 0 production in undisturbed ecosystems. 

Consider a conceptual model of N z O production in a unit area o 
unperturbed land with a given vegetation/soil complex. Litterfall 
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associated with vegetation cycles supplies fresh organic matter to 
the soils. The organic matter is decomposed and mineralized by 
microbial activity, thus delivering N for nitrification and denitrifica- 
tion. The decomposition rate is regulated by soil temperature and 
other parameters of the soil environment. Under most soil condi- 
tions, both nitrifying and denitrifying bacteria are active. The rela- 
tive contribution of N 2 0 produced via the nitrification and denitrifi- 
cation pathways is determined largely by the degree of water 
saturation and aerobicity of the soils. The wetter the soil, the lower 
the oxygen content and the more likely denitrification processes are 
to dominate. Soil water and oxygen status are, in turn, governed by 
the net water supply at the surface, soil water-holding capacity, and 
drainage properties. Other inherent soil properties, such as pH and 
fertility, also modulate the microbial activity and N 2 0 production. 

The controls on N 2 0 production from natural soils can be sum- 
marized as (1) carbon and nitrogen availability, (2) temperature, (3) 
soil moisture. (4) soil oxygen status, (5) soil fertility, and (6) soil pH. 
Because global data on soil pH are of questionable quality, this fac- 
tor was not included in our analysis. 

Global Data for the Toy Model of N 2 0 Production 

We now present a synthesis of our conceptual understanding of 
the controls on N 2 0 fluxes into a quantitative framework, i.e., into a 
model. The ultimate objective of the model is to evaluate what con- 
tribution natural ecosystems make to the global N 2 0 budget and 
how the contribution would change with global change. Given the 
paucity of measurements, our immediate goal is to test the concep- 
tual model of N 2 0 fluxes and identify gaps where measurements 
and analyses are needed. 

The focus of the model requires data or inputs that must be 
global in domain and span at least a year. The primary gridded data 
sets available are described below: 

• Soils; The 1:5M FAO/Unesco soil maps (1971-1982, 1974, 1988) 
include information on dominant soil units, associated soils 
(when they cover at least 20% of the area) and inclusions, texture 
of the dominant soil, and slope. We used information on the 
major soil units (106) and texture of the topsoil (upper 30 cm) 
from the data digitized at 1° x 1° resolution for the globe by Zobler 
(1986). The spatially dominant soil units in the 1° grid cell were 
recorded. 

Texture reflects the relative proportions of clay (particles less than 
2 pm in size), silt (2-50 pm), and sand (50-2000 pm). Three tex- 
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ture classes for the topsoil are distinguished in the FAO soil 
maps: coarse, medium, and fine. At 1° resolution, a spatially 
dominant texture cannot always be established. And so, four 
other texture classes were included in the digital data base: 
coarse/medium, coarse/fine, medium/fine, and coarse/medium/ 
fine. Furthermore, texture information was not available every- 
where in the FAO soil maps, in particular for some of the organic 
soils e.g., Histosols. A separate class, "organic,” was added to the 
texture description. These additional five classes composed <10% 
of the land surface. 

The soil type information is used to translate directly into fertility 
and drainage, while the type and texture information are used to 
derive soil storage capacity. These are described briefly below and 
in detail in Bouwman et al. (in preparation). Figure 1 shows the 
global distribution of soil texture, grouped into six classes. 


Texture 



Figure 1. Distribution of soil texture: 0 = organics, ] = coarse. 2 = coarse /medium, 3 = medium 4 = 
medium/ fine, 5 =fine. Spatial resolution is l°x 1°. 
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• Vegetation: The satellite-derived normalized difference vegetation 
index (NDVI) is a measure of primary productivity of the vegeta- 
tion. The NDVI is the difference between the radiances in the visi- 
ble (0.58-0.68 pm) and near-infrared (0.725-1.1 pm) wave- 
lengths, normalized by the sum of the radiances. The radiances 
are measured by the advanced very high resolution radiometer 
(AVHRR) on board the NOAA series of polar-orbiting satellites 
(Tarpley et al., 1984). In this study, monthly NDVI composites for 
1984 were gridded at 1° x 1° resolution for the globe, and then 
summed to produce the annual total. Although a digital data base 
of vegetation types exists (Matthews, 1983), vegetation informa- 
tion is used only for verifying measurement site characteristics 
and for analysis of model results. 

• Climate: Shea (1986) has produced climatologies of monthly mean 
surface air temperature and precipitation, at 2° x 2° resolution for 
the globe, from the available station observations. In the concep- 
tual model, the climate parameter important for N 2 0 production 
is soil temperature. Lacking a global climatology for soil tempera- 
ture, we used surface air temperatures instead. This would intro- 
duce phase errors, of up to a season in middle and high latitudes, 
in the seasonality of N a O production. 

The finest spatial resolution of these data sets is 1° latitude by 
1° longitude, so that there are -15,000 grid cells for the ice-free 
land surface. The climate data are monthly means. With this reso- 
lution, the model cannot resolve episodic effluxes of N 2 0 after 
rainstorms and spatial “hot spots” which are often reported. The 
importance of such high-frequency, local events in the global bud- 
get has not been established. The model must parameterize, in 
some way, their integrated effect, and evaluate their contribution 
to the global annual flux. 

A Toy Model for N2O Production 

Fluxes of N 2 0 from temperate grasslands have been extensively 
modeled based on the comprehensive measurements there (e.g., 
Parton et al., 1987, 1988). For the other soil types around the globe, 
there is a dire lack of N 2 0 flux measurements, let alone quantitative 
relationships between N 2 0 fluxes and the various controlling para- 
meters. We do not know if the relationships established for grass- 
lands can be extrapolated elsewhere. 

Our strategy for the toy model for N 2 0 production is to first trans- 
late ideas about relative importance of each of the controlling para- 
meters into ranked nondimensional indices. The indices are scaled 
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from 0 to 10, or from 0 to 5, with high numbers signifying impor- 
tance for N 2 0 production and 0 signifying no production. Our choice 
of a scale of 5 or 10 depends on the number of levels we think we 

can discriminate within the scale. 

Using nondimensional variables to denote relative importance is 
not a new concept. It was employed, for example, to represent tem- 
perature and moisture controls in the N 2 0 model of Mosier and Par- 
ton (1985). Such a translation is straightforward for numeric data 
such as temperature. In this toy model, we also created nondimen- 
sional scales for ordinal data such as fertility and drainage by carry- 
ing out a subjective ranking of soil units and texture information. 
The control indices are then combined to form the index for N z O 
fluxes. Comparison of the scaled nondimensional N z O fluxes with 
field measurements provides an evaluation of the validity of the 
model and, one hopes, an algorithm for calculating N 2 0 fluxes from 
the scaled variables. 

The model is illustrated schematically in Figure 2. A central part 
of the problem is to model how soil water regulates the degree of 
nitrification/denitrification and hence N 2 0 production. A bucket 
model of water balance is used. The bucket model takes into 
account differences in soil drainage and topsoil texture and deter- 
mines the water and oxygen status of the topsoil. For both nitrifica- 
tion and denitrification. N 2 0 flux is also proportional to the 
amounts of carbon and nitrogen available in the soil, the rate of 
mineralization of the carbon and nitrogen, and soil fertility. 

With this approach, we need to derive the global distribution of 
the five major regulators of N 2 0 production. Three of them vary 
monthly: decomposition of soil organic matter (SOD), water avail- 
ability (WATER), and oxygen limitation (OXYGEN): the remaining 
two, soil fertility (FERT) and carbon and nitrogen availability (CAR- 
BON), are constant through the year. Below, we describe in detail 
the derivation of each of the regulators and their synthesis to yield 
the nondimensional N 2 0 flux. 

Soil Carbon and Nitrogen Available 

Litterfall and root decay are the major sources of carbon and 
nitrogen to the soil. For most ecosystems, there is an abundance of 
litter on the ground throughout the year. Hence the seasonal varia- 
tion of C and N mobilization in the litter is generally governed by 
decomposition rates rather than by seasonal variations in litterfall. 
Here we assumed that geographic pattern of annual litterfall is the 
same as that of annual net primary productivity of the ecosystems, 
for which the satellite-derived NDVI has been shown to be a good 
correlate (Goward et al., 1986; Box et ah, 1989). Because litterfall is 
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not in phase with productivity, the annual NDVI totals, rather than 

NDVI for individual months, were used. 

Monthly NDVI values range from - 0.1 to 0.5, and the annua 
totals range from -0.1 to 4. To be consistent with the other factors 
used in the study, these NDVI totals were, in turn, rescaled to range 
from 0 to 10 to obtain the index CARBON (Figure 3). We note that 
we could have used, instead of the NDVI, a global distribution of net 
primary production (NPP) obtained by assigning literature values of 
annual NPP to the digital data base of vegetation. The use of the 
NDVI captures the variability of NPP within each vegetation type, at 
the same resolution as the soil data. 

Delivery of Nitrifiable N , v ,, . f 

The rate of carbon and nitrogen delivery is governed by the rate ot 
decomposition and mineralization of soil organic matter, w ic are 
regulated by soil temperature, soil moisture, soil fertility, and soil 
texture. Here, we introduced a factor in the model to represent the 

Carbon + 9{itroyen Availability 
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temperature effect on the supply of nitrogen and carbon through 
decomposition of organic matter (SOD). 

Three temperature dependencies for the factor SOD were investi- 
gated (Figure 4). The first (referred to as SOD1) is an exponential 
function describing the temperature effect on N 2 0 fluxes obtained for 
semiarid grasslands (Mosier and Parton, 1985; Parton et al., 1988). 
SOD1 = 10 at T = 50°C, and SOD1 = 0 for T < 0°C, with a rapid 
increase between 10 and 30°C. The second equation (referred to as 
SOD2) is a quadratic function adapted from Parton et al.’s (1988) 
equation for N 2 0 fluxes from grasslands. The optimum (SOD2 = 10) 
is at T = 35°C. To avoid negative values, SOD2 was set to zero at T < 
7°C and T > 50°C. The third relationship tested (referred to as SOD3) 
is a set of linear functions for four broad ecosystem groups: broad- 
leaved vegetation, needle-leaved vegetation, grasslands, and tropical 
vegetation. The functions were derived from observations of C0 2 evo- 
lution from soils (Fung et al., 1987). They were set to SOD3 = 10 at T 
= 50°C. The different slopes of the four functions reflect the differ- 
ences in litter composition for the broad biome groups. For compara- 
ble temperature ranges, grasslands, having a higher fraction of easily 
decomposed material in its detritus than needle-leaved woody vegeta- 
tion, would have a faster relative decomposition rate. The maximum 



S0D1 " S0D2 • x SOD3 (Needlelvd) 

SOD3 (Broadlvd) — S0D3 (Tr op/subtrop) S0D3 (Grasslands) 


Figure 4. Index of the rates of soil matter decomposition (SOD) as a 
function of temperature. 
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values at T = 50°C do not capture the adaptation of microbes to their 
environmental conditions. 

We used the exponential curve (SOD1) in the reference formula- 
tion for N 2 0 production (see “N 2 0 Production” below). 


Soil Water and Oxygen Status 

A key to distinguishing the pathways of nitrification vs. denitrifica- 
tion is the degree of saturation and aerobicity of the soil, which is, in 
turn, determined by the soil drainage properties, as well as by the 
amount of water present in the topsoil vs. the maximum amount of 
water allowed in the topsoil, i.e., the soil water storage capacity (SSC). 

A first step in modeling the water balance of the topsoil is the 
determination of the soil water storage capacity. SSC is the total 
amount of water held in the upper 30 cm of the soil profile at field 
capacity (soil water potential of 10 kPa), after internal drainage and 
redistribution have ceased. Global distributions of SSC are used in 
general circulation models (GCMs); recent hydrology schemes assign 
SSC based on vegetation and/or soil characteristics. We prescribed 
SSC based on soil texture for most soils. SSC is 120 mm, 80 mm, 
and 40 mm for fine-, medium-, and coarse-textured soils, respec- 
tively. For a number of soils, covering about 30% of the ice-free land 
surface, properties other than texture exert the primary influence on 
water retention; their SSCs were assigned without regard for tex- 
ture. Generally, these are fine-textured Ferralsols and Vertisols 
whose aggregates and swell-shrink properties, respectively, reduce 
their water-holding capacity to that of medium-textured soils. The 
global distribution of SSC is shown in Figure 5. 

Soil Water Budget Model 

The change in soil water content (SWC) is the difference between 
the supply and demand of moisture at the surface. Supply is gov 
emed to a large extent by precipitation, while demand is governed 
by evaporation through soils and transpiration through plants. 
Runoff (on sloping land) or ponding (on level land) occurs when the 
net input of water plus the initial amount of soil water exceed soil 
water storage capacity after drainage has been accounted for. 

There are several soil moisture models used in general circulation 
models, ranging from the simple bucket model of Manabe (1969), 
where SSC was uniformly 15 cm, to the simple biosphere model 
(SiB) of Sellers et al. (1986) and the complex biosphere-atmosphere 
transfer scheme (BATS) of Dickinson (1984) and Dickinson et al. 
(1986), which take into account the effects of biome differences on 
soil water balance. Not only do these recent models distinguish 
among soil evaporation and plant evaporation and transpiration. 
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Soil "Water Storage Capacity 



1 l I 1 1 1 1 . L -J 

O 60 120 ISO 



15 40 60 BO 100 120 


(mm) 


Figure 5. Distribution of SSC. 


they also model explicitly vertical profiles of moisture in the soil. 
These recent models are not yet fully validated in GCMs and are too 
complicated for the toy model at hand. However, in a simple bucket 
model where supply and demand are independent of soil moisture 
itself, the solution to the soil moisture equation is not uniquely 
determined; it is dependent on the initial SWC assumed, unless the 
SWC = 0 or SWC = SSC (runoff) sometime during the year. 

In this study, we chose a soil moisture model that is simple in 
design, and whose solution is not dependent on an arbitrarily cho- 
sen initial condition. We adapted the Mintz and Serafini (1981) 
model for the calculation of the monthly soil water content (SWC m ) . 
In this model, net supply is the difference between monthly precipi- 
tation (P m ) and evaporation from wet canopies (EI m ), while demand 
is the sum of transpiration through plants and evaporation from 
soils (ETS m ) . The time step of a month corresponds to the time reso- 
lution of the global precipitation climatology used. 
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SWC m - SWC m _! + (P m - EI m ) - ETS m 

Potential evaporation (PE m ), a measure of the maximum demand 
for moisture by the atmosphere, was calculated from the surface air 
temperature according to Thomthwaite (1948). Mean monthly sur- 
face air temperatures (T m ) and precipitation (P m ) were obtained from 
the climatologies compiled by Shea (1986). For simplification, pre- 
cipitation in the form of snow was treated as rainfall. Three mois- 
ture regimes were considered, dependent on the relation between P m 
and PE m : 


o 

II 

EI m =0 


ETS m = PE m x p m xa 

< PP'm 

Elm = 


ETS m = (PE m - P m ) x p, 

Pm - P^m 

= PE m 


ETS m =0 


where 


a = 0.4 


and 

Pm = 1 - e-T [{SWC^j + (P m - EI m )/2}/SSC - 5] 

The coefficient a expresses the ratio of the amount of water extracted 
from the topsoil to the extraction from the full rooting zone (as noted 
before, we consider only the topsoils). The function f$ describes the 
maximum water extraction as a function of soil water content and 
soil characteristics. Its parameter y depends on topsoil texture and 
mineralogy, while 8 represents the water unavailable to plants, i.e., 
the intercept of the water extraction curve p. The values of 8 and y 
are dependent on soil texture (see Bouwman et al., in preparation). 
For fine textured soils, y = 6, resulting in a strong decrease in water 
extraction below SWC/SSC « 50%. Due to the selected value of 8 for 
clays, water extraction at SWC/SSC < 10% is reduced. In sands and 
medium-textured soils P sharply declines at values of SWC/SSC of 
about 40% and 20%, respectively. Because SWC is a function of p, 
the monthly equilibrium SWC is achieved independently of the initial 
water content at the start of the simulation. 

Effects of Drainage 

In the calculation of the soil water content of the topsoil described 
above, drainage properties of the soil were not considered. Drainage 
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determines how excess water is removed from the soil, and it is also 
an indicator of the soil aeration. We ranked the soil drainage index 
DRNG based on soil texture and other soil properties, with a high 
rank favoring N z O production. The very poorly drained soils (DRNG = 
5) include those soil groups strongly influenced by groundwater, i.e., 
Gleysols and Histosols, as well as soils with permafrost within 200 
cm of the soil surface (gelic soil units). Freely drained soils, such as 
Xerosols and Yermosols, occurring in deserts, were assigned a 
drainage rank of 1 . The distribution Is shown In Figure 6. 

We are uncertain how to model quantitatively the effects of 
drainage on soil water content. Adding diffusion at the base of the 
“bucket" may be an appropriate approach, but the values of the dif- 
fusion coefficients for different drainages are not well known. Here, 
we introduced the factor soil water status (SWS), an index of the soil 
water saturation when drainage effects are considered (Table 1). It is 
clear that the highest SWS rank of 10 would be assigned to a poorly 
drained soil when the monthly soil water content SWC is close to 
the storage capacity SSC of the soil. It is not clear how to rank the 
intermediate values of drainage and SWC/SSC. We noted several 
points: (1) Distinguishing saturation levels <20% is not important, 
so that the SWC/SSC scale can be divided into ten levels as shown 
(2) We linearly increased the SWS scale up the saturation scale. (3) 
We filled in the rest of the table by assuming that N 2 0 production 
would likely be asymptotic at high saturation and poor drainage. 
There is thus much arbitrariness in the SWS scale. It nevertheless 
represents a first attempt at quantifying our descriptive understand- 
ing of the effects of soil drainage characteristics on soil moisture. 

Nitrification and Denitrification Potentials 

It remains to describe the effects of soil water status on oxygen 
suppiy in the soils, i.e., on nitrification and denitrification. Two fac- 
tors, water availability (WATER) and oxygen limitation (OXYGEN), 
were derived as indices for nitrification and denitrification poten- 
tials, respectively. These two indices are based on the two moisture 
indices derived by Parton et al. (1988) for the two processes. 

The influence of soil water status on nitrification and denitrifica- 
tion has been documented by many authors (Mosler and Parton, 
1985; Klemedtsson et al., 1988; Groffman and Tiedje, 1988). Aerobic 
microbial activity Increases with soil water content until a point is 
reached where water displaces air and restricts oxygen diffusion. 
Maximum rates of nitrification occur at the highest water content at 
which soil aeration remains nonlimiting, with SWC -60-80% of SSC. 
And so we assigned maximum nitrification potential WATER at SWS 
around 6-8 (Table 2). 
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Figure 6. Dtelrtbuflon of the drainage Index DRNC. ’ Value s of I and 5 denoK pood- and poor 
drainage soils, respectively. 


Table 1: Soil water status , as a Junction of the soil drainage scale and calcu- 
lated soil water content/ soil water storage capacity 


swc/ssc 

(%) 


Drainage Scale 


0-20 
20-30 
30-40 
40-50 
50 - 60 
60 - 70 
70 - 80 
80 - 90 
90-100 
> 100 % 
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Table 2: Water availability as a function of soil water status in current and preceding months 


SWS, 

preceding 

month 

1 

2 

3 

SWS, current month 
4 5 6 7 

8 

9 

10 

1 

1 

2 

4 

6 

9 

10 

10 

10 

6 

1 

2 

1 

1 

3 

5 

8 

9 

10 

10 

6 

1 

3 

1 

1 

2 

4 

7 

8 

9 

10 

5 

1 

4 

1 

1 

2 

4 

7 

8 

8 

8 

4 

1 

5 

1 

1 

2 

4 

6 

7 

7 

7 

3 

1 

6 

1 

1 

2 

4 

6 

7 

7 

6 

2 

1 

7 

1 

1 

2 

4 

6 

7 

7 

6 

2 

1 

8 

1 

1 

2 

4 

5 

6 

6 

5 

2 

1 

9 

1 

1 

2 

3 

4 

6 

6 

4 

2 

1 

10 

1 

1 

2 

3 

4 

5 

6 

4 

2 

1 

No dependence' 

* 

1 

2 

3 

5 

7 

10 

10 

5 

21 

*Scalars depend only on the SWS in the current month. 


The maximum N 2 0 production by nitrifier denitrification is 
between 80 and 100% water-holding capacity, while the peak for 
respiratory denitrifier N a O production may be at water contents 
exceeding 100% water-holding capacity (Klemedtsson et al„ 1988). 
Hence OXYGEN, the index of oxygen limitation and denitrification 
potential, was assigned a maximum value for SWS = 10 (Table 3). 

Field measurements generally show high pulses of N 2 0 efflux from 
the soils after rainfall. Lacking any statistics, we do not know how 
significant these episodic pulses are in the global N 2 0 budget. Fur- 
thermore, there are presently no gridded climatologies of daily rain- 
fall even for a continent. To test the Importance of wetting and diying 
using the data available, we made the nitrification and denitrification 
indices dependent on the soil water status of the previous month in 
addition to the current month (see Tables 2 and 3). We hypothesized 
that dry soils that are wetted (large difference in SWS between previ- 
ous and current months) are more favorable for N 2 0 production than 
soils with constant SWS. In this way. pulses of N 2 0 at the onset of a 
wet season would give a higher average monthly N 2 0 flux than in wet 
months preceded by moist conditions. The hysteresis effect observed 
(Groffman and Tiedje, 1988) is Included by making the Indices under 
drying conditions lower than those under equal amounts of wetting. 
This procedure imposes a one-month memory on the system, and the 
memory is probably much longer than that represented in the 
episodic pulses observed. This is thus one of the most uncertain 
aspects of the toy model, but an aspect that can be tested with long 
time-series measurements such as that of Parton et al. ( 1 988) . 
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Table 3: Oxygen limitation as a function of soil water status in current and preceding months 


sws, 

preceding 

month 

1 

2 

3 

SWS, current month 
4 5 6 7 

8 

9 

10 


\ 

1 

1 

1 

4 

6 

8 

10 

10 

10 

10 


2 

1 

1 

1 

4 

6 

8 

10 

10 

10 

10 


3 

1 

1 

1 

3 

5 

7 

9 

10 

10 

10 


4 

1 

1 

1 

3 

4 

6 

8 

9 

10 

10 


5 

1 

1 

1 

2 

3 

5 

7 

8 

9 

10 


Q 

1 

1 

1 

1 

2 

4 

6 

7 

8 

9 


7 

1 

1 

1 

1 

2 

3 

5 

6 

7 

8 


Q 

1 

1 

1 

1 

2 

3 

4 

5 

6 

7 


9 

1 

1 

1 

1 

2 

3 

4 

5 

6 

7 


10 

1 

1 

1 

1 

2 

3 

4 

5 

6 

7 


No dependence' 

* 1 

2 

2 

3 

3 

4 

5 

6 

8 

10 


•Scalars depend only on the SWS in the current month. 


In the reference case for N 2 0 production (see “N 2 0 Production”), 
we included the WATER and OXYGEN dependence on the water sta- 
tus of the previous month. 


Soil Fertility 

The last control on N 2 0 production we need to quantify Is soil fer- 
tility. Matson and Vitousek (1987, 1990) have shown a strong rela- 
tionship between N 2 0 flux and soil fertility in the tropics and sub- 
tropics, and that relationship is useful for regional estimation of 
N 2 0 emission. In this study, soil fertility for each of the 106 soil 
units was ranked subjectively based on general understanding of 
the cation exchange capacity and other soil properties. Because we 
did not feel that we were capable of distinguishing as many as ten 
levels of soil fertility, five ranks were used. A high value denotes a 
fertile soil, conducive to high levels of N 2 0 production. For example, 
Ferralsols, occurring mainly in the tropics, are strongly leached soils 
and their fertility is low because of their low cation exchange capac- 
ity, presence of alumina on the exchange complex, low content of 
weatherable minerals, and phosphorus fixation. A rank of 1-2 was 
assigned. In contrast. Chernozems and Kastanozems, underlying 
much of the agricultural lands in midlatitudes, are fertile soils to 
which a rank of 5 is assigned. The global distribution of the index 
FERT is shown in Figure 7. 

N 2 0 Production 

With the controls on N 2 0 fluxes quantified, it remains to combine 
them to yield PROD, the nondimensional index of N 2 0 production. 
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Figure 7. Distribution of the fertility index FERT. High values denote fertile soils , conducive to N 2 0 
production. 


We note that the five factors chosen are not independent. In particu- 
lar, CARBON, scaled from the annual NDVT integral, captures geo- 
graphic variations in temperature and soil moisture, as well as vari- 
ations in soil fertility. While it may be argued that annual N 2 0 
production may be proportional to CARBON alone, we still need to 
obtain seasonal variations in the production. 

There are many ways to combine the factors. Lacking any global 
information, we first assumed that all five controlling factors are of 
equal importance, i.e., the maximum fertility factor has the same 
effect as the maximum oxygen limitation factor as far as N 2 0 produc- 
tion is concerned. Hence, although FERT was scaled from 1 to 5 
because of our inability to discriminate further, FERT was multiplied 
by two to normalize to the other factors. The effects of the five factors 
may not be equal and can be analyzed in sensitivity analyses. 

We then modeled the nondimensional N 2 0 production every 
month as the geometric mean of all five controlling factors. We 
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chose a geometric mean rather than an arithmetic mean because a 
factor of zero, such as when temperature falls below 0°C, automati- 
cally shuts down N 2 0 production. Also, a low value for one of the 
factors would lower the N 2 0 production index. For example, values 
of 1 and 9 for two factors are given less weight than 5 and 5, which 
yield the same arithmetic mean. 

The degrees of nitrification and denitrification are both governed 
by the water and oxygen supply in the soil, and the products of nitri- 
fication provide the substrate for denitrification. However, the 
processes themselves are carried out by different microbial popula- 
tions and are, to a large extent, independent of each other. It is thus 
not correct to multiply the denitrification and nitrification indices 
OXYGEN and WATER. In their model of N 2 0 production in the Colo- 
rado grasslands amended with urine, Parton et al. (1988) summed 
the contribution of the two processes to obtain the total N 2 0 produc- 
tion. There, sufficient data allowed the relative contributions of the 
processes to be obtained as a constant in the regression equation 
between the model and the flux measurements. Parton et al. found 
that nitrification contributed 60-80% of the total flux. This is in 
sharp contrast to tropical soils where denitrification dominates (e.g., 
Matson and Vitousek, 1990). We do not have information on the rela- 
tive importance of nitrification and denitrification for all ecosystems, 
nor do we have sufficient N 2 0 flux measurements to carry out a 
regression analysis similar to that of Parton et al. (1988). Thus, this 
version of the toy model cannot discriminate between N 2 0 produced 
via the different processes, and the product of OXYGEN and WATER 
can only be viewed as an index of the total production. 

The monthly N 2 0 production index, PROD, is calculated as fol- 
lows in the reference case: 

PROD = [OXYGEN x WATER x SOD x FERT * xCARBONl 1/5 

where PROD, OXYGEN, WATER and SOD are indices calculated 
monthly: FERT* = 2 x FERT; and CARBON has a fixed value for all 
months. 

Model Evaluation 

The model was applied at each of the 1° x 1° grid boxes for the 
globe, to yield, every month, the global variation of N 2 0 production. A 
series of sensitivity experiments was carried out. Instead of the refer- 
ence exponential temperature dependence (SOD1), the quadratic 
(SOD2) and linear (SOD3) curves were used separately. Other sensi- 
tivity experiments selectively excluded each of the controlling factors. 

Our next step is to compare our modeled indices for N 2 0 produc- 
tion from the reference case and each of the sensitivity experiments 
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with fluxes directly measured in the field. With the comparison, we 
hope to establish that the geographic and seasonal variations in the 
nondimensional N 2 0 indices capture those found in the field. We 
carried out a literature survey of measurements of N 2 0 fluxes from 
natural systems. Only six different ecosystems (wetlands, temperate 
forests, steppes, tropical savannas, tropical dry forests, and tropical 
rain forests) are represented. Some of the measurements spanned 
less than a month. 

Where measurement conditions corresponded with the vegeta- 
tion/soil/climate conditions in the digital data bases used as inputs 
to our toy model, we compared the measured N 2 0 flux and the mod- 
eled PROD, both averaged over the period of the measurements. For 
measurement sites where discrepancies exist between field condi- 
tions and model inputs, we recomputed PROD using the conditions 
of the measurement sites. 

By averaging the data and model results over the period of mea- 
surements, we hope to ameliorate the impact of the phase errors 
associated with using surface air rather than soil temperatures in 
our analysis. Small-scale variability may still confound the compari- 
son between point measurements and model calculations at 1° reso- 
lution. Ideally, we would average measurements in the 1° grid box 
for comparison. This being impossible, we averaged the modeled 
indices over a few grid boxes adjacent to the measurement sites, in 
hopes that variability within a grid box is reflected in variability 

among several grid boxes. A total of 30 data points resulted for the 
comparison. 

Of all the sensitivity experiments, the standard model as pre- 
sented above yielded the highest correlation coefficient This is 
shown in Figure 8 together with the best line fit to the data A qua- 
dratic fit 4 


y = 78.9-52.3x + 12.4x 2 

was found to produce the highest r 2 = 0.57 between the model and 
observations. Here x represents the modeled monthly averaged 
index PROD, and y the monthly averaged N a O flux measured (g 
N/ha/mo). 6 

Figure 9 shows the global distribution of PROD, summed over a 
year. The highest N 2 0 production is found in the equatorial regions 
where both temperature and NDVI exhibit their maximum values. In 
the reference case, the equatorial regions (30°N-30°S) account for 
80% of the global production. While midlatitude Chernozems and 
Kastanozems have the highest fertility, and permafrost regions have 
highest oxygen limitation, their inclusion only reduces the latitudi- 
nal gradient in N a O production. 


„ M ItiUltciU UUUUOii i 4-XiHlu., 


258 


Modeling the ‘Earth System 



Figure 8. Comparison of PROD, the modeled N 2 0 index, with measured 
N 2 0 fluxes (boxes). The thick and dotted lines represent the best line fit ± 1 
standard deviation. 


Conclusion 

The regression equation provides the transformation of the mod 
eled nondimensional N z O indices into dimensional N 2 0 fluxes. The 
global N 2 0 emission from natural soils obtained is 7 Tg N/yr (wit a 
confidence interval of 3.1 to 13.4 Tg/yr). That our modeled estimate 
is within previous estimates for this source should be no surprise, 
as all extrapolation schemes are based on the same limited set of 
measurements of N 2 0 fluxes. Of the -15,000 1° x 1° ice-free land 
cells only 30 are represented with measurements. Thus the toy 
model, while more elaborate, does not provide a more accurate esti- 
mate of the global N 2 0 emission than other extrapolation schemes. 
However, the model and the sensitivity experiments do present a 
framework for analyzing how different controls affect N 2 0 produc- 
tion in different regions and for identifying the gaps in our under- 
standing of N 2 0 emissions from natural soils, as well as a strategy 
for measurements. 

A major implication of the toy model is that the tropics account 
for >80% of the global emission, within the range inferred by Prinn 
et al (1990) from the atmospheric N 2 0 distribution. The large tropi- 
cal source, in our model, is a result of the equal weighting given to 
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Figure 9. Distribution oj nondimenslonal index oJN 2 0 productior 


80 

summed over the year. 


1 hC latitudinaI gradient in N 2 0 production is 
dominated by the controlling factors CARBON (carbon and nitrogen 
availability, proportional to NDVI) and SOD (decomposition rates 

r r ° P nul^riL t0 temperatUreJ ' 7116 two moisture indices for nitrifica- 
tion (WATER) and denitrification (OXYGEN) show no distinctive lati- 
tudinal gradients, while soil fertility peaks in the Chernozems and 
Kastenozems of midlatitudes. Whether the latitudinal gradient of 

2 production is as steep as that estimated by the toy model can 
e assessed by a detailed comparison and analysis of flux measure- 

merits from different locations. 

We note Matson and Vitousek (1990) appear to have a significantly 
lower emission from the tropics than that estimated by the toy 
model. The difference may not be real, but may lie in the choices of 
areas or soil types included in the two calculations or in land use 
effects not included in this study. With the series of geographic data 
bases used m the toy model, it is straightforward to cany out parallel 
analyses so that the discrepancies can be evaluated and understood. 
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Improvement of the global estimate of N 2 0 emission from natural 
soils clearly requires more flux measurements. The toy model sug 
gests that measurements should span at least a complete annual 
cycle and should be carried out along several transects that cover 
di^ent gradients in temperature, precipitation, soil texture and 
fertility and vegetation productivity. It is also clear that complete 
life description with soil, vegetation, and climatic data, is cruoal 
for understanding the controls on N z O production. 
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